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Abstract 

Constant Proportion Portfolio Insurance (CPPI) is a strategy designed to give participation in 

a risky asset while protecting the invested capital. Some gap risk due to extreme events is often 
kept by the issuer of the product: a put option on the CPPI strategy is included in the product. In 
this paper we present a new method for the pricing of CPPIs and options on CPPIs, which is much 
faster and more accurate than the usual Montc-Carlo method. Provided the underlying follows a 
homogeneous process, the path-dependent CPPI strategy is reformulated into a Markov process in 
one variable, which allows to use efficient linear algebra techniques. Tail events, which are crucial 
in the pricing are handled smoothly. We incorporate in this framework linear thresholds, profit 
lock-in, performance coupons... The American exercise of open-ended CPPIs is handled naturally 
through backward propagation. Finally we use our pricing scheme to study the influence of various 
features on the gap risk of CPPI strategies. 
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1 Introduction 



Constant Proportion Portfolio Insurance is a dynamic strategy designed to give participation in some 
risky asset while protecting the invested capital f Mer71] . This is achieved by rebalancing periodically 
between a risk-free asset (Zero-Coupon bond) and a risky asset(share, index, fund, fund of funds. . . ). 
In the simplest form, if the underlying asset has no jumps and if one can rebalance continuously, the 
final payoff depends in a deterministic way on the risky underlying. However both hypothesis are 
strong and do not fit real market conditions. Relaxing them, the strategy is not as efficient: there is 
a small chance of not recovering the whole capital invested at maturity. This gap risk may be hold by 
the issuer, so that the principal is really guaranteed to the investor. In this case, there is an option 
included in the product, which must be priced and hedged. With a discrete rebalancing scheme, there 
is a closed formula for this option price if the underlying follows a Black-Scholes diffusion |BBM05] . 
Moreover, this formula can be generalized to jump-diffusion models, and more generally to Levy 
processes. However, this works only for an idealized CPPI product where there are no spreads on 
the risk- free and financing rates, no fees, a natural bond floor... Unfortunately, real CPPIs always 
have such features, which prevent from using the closed formula. As CPPIs are very path dependent, 
they are usually priced through Monte-Carlo simulations (see [BK95J for an example). Extreme value 
theory has also been used to estimate the gap risk of such products jBP02j . 

Monte-Carlo pricing is perfectly suited to path dependence but the dimension of the problem is 
generally quite large. For example, a monthly CPPI defined on a single risky asset with monthly 
rebalancing and a maturity of 20 years requires 240 fixings. This means a 240-dimensional Monte- 
Carlo integration. Furthermore, the tails of the distributions play important roles in the pricing. 
First, the lower tail gives the gap risk part. It must be computed with enough precision to produce 
a reasonable Put price and an accurate Delta. The upper tail is important to reach the correct mean 
value of the CPPI strategy. The payoff distribution is close to a shifted lognormal distribution with 
high leverage: most trajectories will end far below the mean and will be compensated by few very high 
terminal values. Control variate can solve part of the convergence problems for the strategy mean 
value, by comparing the payoff of the real CPPI to the payoff of a self-financing CPPI. However that 
does not solve all convergence issues: in particular the optional part, the gap risk, is not smoothed by 
basic control variate. The high path dependence and barrier-like structure of the strategy makes the 
Put price and its Delta converge very slowly. Moreover, more and more CPPIs are open-ended: the 
investor can exit at any time. The Monte-Carlo method is not well suited to this American feature. 

We propose here an alternative way for pricing CPPI strategies and derivatives. Instead of sim- 
ulating trajectories, transition probabilities are computed. This takes care smoothly of distribution 
tails and produces smooth greeks. Initially, there are four correlated variables to follow: the risky 
asset, the CPPI strategy value, the risky asset exposure and the guaranteed amount. In many situa- 
tions, the problem can be reformulated in such a way that it becomes one-dimensional and Markovian. 
Instead of diffusing the underlying, the CPPI index value is evolved. Transition matrices between 
two rebalancing dates are computed. From the Markov property, the final distribution can be easily 
computed from the starting point and one can even compute the full transition matrix. The main 
advantages of this technique are its speed and the smoothness of greeks. For example in the simplest 
cases, only a few milliseconds are needed to get an accurate price. More complex cases with profit 
lock-in, artificial cushion or coupons take more time but never more than a few seconds with cor- 
rect approximation. As the pricing is done through backward propagation, products with American 
optionnality such as open-ended CPPIs are priced naturally. 

After setting the model assumptions in section [2] our pricing method is described in section |3j 
Features which deserve a special treatment are considered in |4[ performance coupons, profit lock- 
in, conditional rebalancing. . . We discuss the influence of the discretization scheme and how to keep 
under control the numerical artifacts arising from the discretization process in section [5} Numerical 
results are presented in the two last sections. Convergence results are shown and compared to a 
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Monte-Carlo pricing in section |6j Then we illustrate the efficiency of the pricing by studying the 
impact of various CPPI features on the final distribution and especially on the gap risk in section 
[7j Finally, two appendices give useful formulas for the jump-diffusion model we use and for vanilla 
CPPI prices. 

2 Assumptions 

2.1 Numeraire choice 

We consider two cases. 

1. Deterministic interest rates. The numeraire is the monetary unit. All features described in this 
paper can be handled in this case. 

2. Stochastic interest rates. In this case we choose the forward risk- neutral measure: the numeraire 
is a Zero-Coupon bond of maturity larger than the CPPI maturity. The assets which are diffused 
are of the form ^^^^ j,^ . The stochastic nature of interest rates can be forgotten provided every 
quantity can be expressed simply in term of this numeraire. Linear threshold and fixed rate 
threshold are therefore excluded, because they have "absolute" value in monetary unit and 
are not multiple of a Zero-Coupon bond. With this numeraire, the threshold would become 
stochastic. It is not possible to reduce the diffusion to a single asset, one has really to model 
two assets. To be able to use stochastic interest rates in this framework, the threshold must be 
defined as a Zero-Coupon curve, with the possibility of adding some spread over the risk-free 
rate. 

Formally, we will write deterministic interest rates, which may come from stochastic rates in 
forward risk-neutral measure. Just remind that this latter case is not compatible with linear or fixed 
rate threshold. 

2.2 Underlying process 

In order to be able to model directly the CPPI strategy value without following the underlying asset 
in addition, we make the assumption that the risky underlying evolution is homogeneous: AS/S does 
not depend on S. However it can depend on time-dependent parameters. 

Black-Scholes diffusion will work. Similarly, models with jumps which respect the Levy property 
(independent increments) fit naturally in the framework. In the appendix, we describe the Kou model 
with exponentially distributed positive and negative jumps but other jump distributions can be used 
(Merton's for instance). The choice of the model depends on the context: for pricing and hedging 
a market calibrated model has to be used, whereas for portfolio optimization or risk management 
historical probability can be used |CT07| . 

The hypothesis on the underlying evolution mainly excludes local volatility models. Stochastic 
volatility models can be incorporated in this framework but the volatility must be included in the 
variables space, so a few volatility regimes would be more appropriate numerically than a continuous 
volatility diffusion |Alb07] . If the CPPI value is discretized on a grid with N points, adding 
points for volatility will replace N x N transition matrices with NN^ x NNa- matrices. Starting from 
N ~ 200, a low N(y, between 3 and 5, will keep the matrices tractable: NN^ ~ 1000. 

3 CPPI Pricing 

In a direct approach, four variables have to be modeled: the risky asset value, the CPPI index value, 
the risky asset weighting and the guaranteed amount. In most cases the problem can be expressed 
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as a Markov chain on one single process, the CPPI index value rescaled by the threshold value. 
Transition operators between two rebalancing dates can be constructed. Their multiplication gives 
the transition operator over the whole CPPI life which is used to price the CPPI strategy or options 
on the strategy. 

In this section we deal with the basic case where there is no profit lock-in, coupons or American 
exercise. These features are addressed in the following section. 



3.1 Definitions 

3.1.1 Self-financing CPPI 

The CPPI strategy starts at to = by investing a nominal amount that we take equal to 1 to simplify 
formulas and ends at t„ = T > 0. The portfolio is rebalanced at times ti for i = 0, . . . , n — 1 such 
that to < ti < . . . < tn- The value of the portfolio at time U is denoted d = C{ti). Similarly, we 
write Si for the price of the risky asset on that date and Fi for its forward price with maturity T. Bi 
denotes the value of the risk-free asset: for a self-financing CPPI, Bi is given by the money market: 

5iocZC(ti,r) (1) 

where ZC(-,T) stands for the Zero-Coupon of maturity T. 

In a first step, the guarante^ Gi is constant and such that Gi = G = 1 for all i. The threshold 
Hi is a deterministic function of the guaranteed amount Gt and of the time U . The natural choice is 

Hi = GZCiU,T) 

The rebalancing criterion at ti states that a stochastic proportion Wi = w{Ci, Hi) of the portfolio 
value should be invested in the risky asset S while a proportion 1 — Wi is invested in the risk-free 
asset B. The proportion Wi is called the risky asset weighting and is usually chosen of the form 

W, = w{G,, m)=m( ^ = m( ^^^) ^ (2) 



Ci J \ Ci / Hj 

where m > is the multiplier of the CPPI. We could also include minimum and/or maximum 
exposures, a cushion limit which triggers full investment in the riskfree asset or any other feature, 
even time dependent, the only restriction being that Wi should be a non negative deterministic 
function Wi = w(A'j) of the ratio 

Ci 

More advanced strategies could be handled at the cost of also following the variable Wi which makes 



our problem 2 dimensional (see section 4.3) 



Prom now on, we refer to the CPPI we have just described as a self-financing CPPI. Closed 
analytical formulas for this special case are derived in appendix [B| 



3.1.2 Spreads, artificial thresholds, fees. . . 

We introduce the following possibly time-dependent spreads, expressed in continuous mode|^ 
• A spread s+ on the risk- free asset, it is usually non-positive. 

^Note that whereas the previously defined processes were in fact continuous, the ones we now introduce are only 
relevant at times ti. 

^If rates and spreads are defined in some other compounding mode, as linear mode for example, the continuous 
spread will slightly depend on the risk-free interest rate r. 
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• A spread s_ on the financing component (i.e. when the risky asset weighting is greater than 1). 
Usually non-negative. 

This means that over the period from tj to tj+i, we have 

where s± = s+lwi<i + s-lwi>i- 

We also introduce a (possibly time-dependent) spread sh for the threshold rate, such that 

If Sh is not zero, the variable Xi = Ci/Hi is no longer the forward value of the CPPI. As the 
threshold is only defined at rebalancing dates, a time dependent spread can be used to "simulate" a 
linear threshold. As for the risky asset weighting, the sole restriction on Hi is that is a positive time 
dependent deterministic function of the guarantee such that the threshold at maturity Hx equals to 
the guarantee. 

Finally, we introduce some managing fees depending on the index value at the beginning of the 
periocQ the risky asset weighting. . . We consider here several types of fees (possibly time dependent): 

• proportional fees fpCi 

• proportional fees replaced with some discounted defeasance rate fd when all the investment is 
invested in the risk-free asset: {fp^Wi>o + fd^Wi=o) Ci 

• fees on the risky investment frWiCi 

• fixed fees // 

The sum of all fees is written as /(Cj, Wi)Ci for convenience. 



3.2 Transition matrix over one period 

The pricing is performed by computing the transition matrix of the CPPI value between the compu- 
tation date to its maturity. For this operator to be computed, the first step is the computation of 
transition matrices over every period. 



3.2.1 Self-financing CPPI 

We consider here the simplest case where there are neither spreads nor fees {s± = sh = f = 0) and 
no profit lock-in {Gi = G = 1). Consider the "rescaled" CPPI value Xi = Gi/Hi. It is a Markov 
process in discrete time and we construct a local (e.g. from one rebalancing date ti to the following 
one) transition operator Mi as follows. The CPPI variation between ti and tj+i is given by 

a+i = w{Xi) Q %i + (l - w{X^))g, ^ 

Oi \ / Di 

In terms of the variables Xi and denoting by Fi the value of the forward on the risky asset with 
maturity T, this looks even simpler: 

^ = 1 + ^X0 (^-l) (3) 



■^One could do the same with fees depending on the value at the end of the period. 
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The evolution operator for the period is defined as 

Mi{x,y)=¥.[5{Xi+i-y) \Xi = x\ 

i.e. it is the density of Xj+i conditionahy on Xi = x. We can compute M.i through the previous 
equation as 

Fi+l 



Mi{x, y) = — P [Xi+i <y\Xi = x] = 



< L{x,y) 



with 



L{x,y) = 1 + 



y-x 

'w{x) x 



3.2.2 General case 

The only difference with the previous case is that the computation are a bit nastier. The value of the 
CPPI portfolio at time tj+i after the fees were paid is 

^ = w{X,)^ + (1 - w{X,)) %i - fia, w{X,)) 



In term of the variable X and the forward F, this reads 

Xi+l Cj+i Hi _ su At ( ^ ,1 V \ -^i+l 



Xi Ci -ffj+i 



w{x,y-^ + (1 - w{x,)) j^r^^-^ - e-^^^^'^^^'mM^,)) 



We invert this relation to get the variation of F conditionally on Xi needed to achieve some value of 
Xi+i/Xi as 



F 



1 



" [1 - w{Xi)) e-"^ * +e f{Ci,w{Xi)) 



e-"-i 



F, w{Xi) y X 

As above, we get the transition kernel between times ti and tj+i as 



Mi{x,y) = |-P 
dy 



Fi+i 
~F~ 



< L{y,x) 



where 



L(x, 



y) = — -— e-"' (1 — w[x))e-'^^ 

w[x) \ X ' 



*+ls±dt , - P+^rdt 

_l_ g Jt,- 



^{x.wix)) 



In the numerical implementation, we make use of the expected value of XiJ^\ conditionally on Xi 
and Xj+i < y. We quote here the result 



IE [Xj_|_ilx,4.i<j/ I Xi 



x\ 



- SH dt , 



(x)P2+ ( {\-w{x))J'^^ -e--^'^^''''^'f{x,w{x)) )Pi 



which is expressed in term of the cumulative Pi 



F 

< L{x,y) and the "conditional" expected 



value P2 = E 



<L{x,y) 



of (see appendix 



A). 
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3.3 Transition matrix over the CPPI life 

The evolution kernel over the whole CPPI life is defined as 

/C(x,y)=E[5(X„-y)|Xo = x] 

allows to compute the expected value of a payoff P{Cn) (as there is no profit lock-in, here we have 
Gi = G = 1 and thus X„ = Gn) as 

E [P{Xn) I Xo] = j dy /C (Xo, y) P{y) 

In operator terms, this is the result of the action of the operator /C on the function P, the result 
being evaluated at Xq: 

E [P{Xn) I Xo] = (/C • P) (Xo) 

Finally we define the partial evolution operators between times ti and tj as 

IC^jix.y) =¥.[6{X, - y)\Xi = x] 

The CPPI process is Markovian: the evolution operator over all periods is given by the product 
of evolution operators over single periods: 

Kij = MiMi+i . ..Mj-i 

where A4„ is the transition operator over the period we considered in the previous section. 

The total evolution is given by the operator 

lC = lC^n=MGMl...Mn-l (4) 

If all local evolution operators Mi are equal to a common operator 7W, this reduces to the power 

/c = /Co„ = Ae 

The fair value today of the product is the expected value 

(/CP)(Xo) = {MoMi . . . Mn-iP) (Xo) 

3.4 Greeks 

Before the CPPI starts, and hence before the first rebalancing of the portfolio, there is of course no 
sensitivity to the risky asset price. Gamma and delta will be identically null until the strategy really 
starts. At a time t such that to < t < ti, consider the spot St and the forward Ft of the underlying 
asset. We can decompose equation ^ as 

and one gets that transition matrix for the remaining time of the first period is given by Mq{x, y) = 
with 



d_ 

dy 



^ < Lo{x,y) 



T f \ ^0 f y-x 

Lo{x,y) = — — — + 1 
Ft \w[x) X 

One can then compute the derivative of the kernel M.o{x,y) with respect to Ft (or the spot St) and 
get the sensitivity of the price to the risky asset as 

AE[P,X.,|X„]./d.5^Pfa) 

with, using Q, 

d}C{Xo,y) _dMo 

dFt -^-^i----^-! 

as only the first kernel depends on the spot of the risky underlying. 
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3.5 Gap Risk 



As sold CPPI are usually guaranteed, we define three useful measures of risk in that context: the 
put price (with a strike equal to the guaranteed amount at maturity), the gap proportion and the 
expected gap loss. The gap proportion is defined as the probability that the final value of the CPPI 
will be below the threshold and the expected gap loss as the expected loss knowing that the CPPI 
ends up below the threshold. The gap proportion is the price of the undiscounted digital put and 
expected gap loss is the quotient of the undiscounted put price by the gap proportion. 



4 Additional features 

CPPI products have often additional features: coupons, profit lock-in, rebalancing conditions... We 
show in this section how they can be handled in the pricing process. 

4.1 Coupons 
4.1.1 Definitions 

The investor might want to be remunerated for the performance but not wait till the CPPI maturity. 
At times t^^i), the investor will then receive money in the form of coupons. To simplify notation, 
we use capital indices to designate coupon dates: Xj means -^c(/) where c(/) G {0, . . . , n} and 
I € {0, . . . , ric} with ric the number of coupons. Two types of coupons have been devised: 

1. Fixed coupon. Fixed coupons are mostly (time dependent) fixed fees except for the fact that 
they are paid to the investor and hence negative. 

2. Performxince coupon. Periodically some proportion q of the performance over the last period is 
paid to the investor. If CJ is the value of the CPPI at time tj immediately before the coupon 
is detached, the value Kj of the coupon paid at for the period (I — 1, /) is computed as 

Kj = q {Cj — C:-i)^ lvKj_i>o 

and 

Ci = CJ - Ki 
Setting XJ = Cj /Hj, the previous equation reads 



Xi = Xj — q ^Xj JiJ^'^^~^ lvK/_i>o 

If no coupon is paid Xj = Xj . In the case a coupon is paid, Kj > and 



l-q\ Hi 



It is also useful to compute the value of the coupon in proportion of the guarantee as a function 
of XJ and Xi: 

— = qlXje -Xi^ie 1 lwj_,>o 



9 



4.1.2 Pricing with performance coupons 

Fixed coupons do not exhibit particular difficulties as they are essentially (negative) fixed fees, hence 
we concentrate here on the pricing of performance coupons. 

Consider A4/_i = K-i-ij be the transition operator, without any coupon, between tc{/-i) ^-iid 
tc(i)- The expected value of the coupon capitalized up to maturity T, conditionally on is given 

by 



ci{x) 



ZC{tj,T) 

Q^w{x)>0 



E 



G 



Xi-i = X 



dz z) I ze 



ZC{tj,T) 

In order to account for the payment of the coupon, we define a rescaled operator Ai: 



(^z-xe^^i^''+'"'''y^ Mi-i{x,z) 



Mi-i{x, y) = J dzS \^z-y - gl^„(^)>o [ z - xe 

By construction this operator satisfies 

Mi-i{x, y)=E [5{Xi - y) \ Xi^i = x] 

which enables us to write /C = K.qNc for the pricing kernel, with 

ICij = Mi Mi+i ... Mj-i 

However, to get the correct price for the strategy, the value of the coupons to be paid must be 
included. The total value of all coupons is 



ric ^ Uc ^ 



G 



gJ^(/Co,/_i-c,)(Xo) 



1=1 



In fact, this hides a semi-direct product structure which, numerically speaking, will be handled 
as matrix multiplication in section [5] For a pair constituted of an operator O and a function /, the 
action of the pair on a function g is given by 

{x) = {0-g) (x) + /(x) 

The (semi-direct) product of two pairs is 

(Ol,/i) (02,/2) = (Ol02, Oi-h + h) 

We claim that the price of strategy with payoff of the form GtP{Xt) including the coupons is given 

by _ _ _ 

G (A^O, Co) {Mi,Ci) . . . {Mn^,Cn}j ■ P 

taken at Xq with the convention co(x) = 0. 



10 



4.2 Profit lock-in 



4.2.1 Definitions 

In order to lock some part of the profit already made, the guarantee G is turned into a stochastic 
process Gj. 

We denote by a minus exponent {XJ) the value of a variable before applying the lock- in. We 
parametrize the lock-in by a multiplicative factor / which relates Xjj^i to XJ_^-^ through 



^i+i _ Gi+i 
Xi+i Gi 

There are two main kinds of profit lock-in: 



fi+i{Xi,X^,-^) 



1. Periodic lock-in. Periodically (every year for example) some fixed proportion p of the perfor- 
mance since the last lock-in is added to the guarantee. As in the case of coupons, we use capital 
indices to designate lock-in dates: Xj means where /(/) G {0, . . . , n} and / G {0, . . . , ni} 
with ni the number of lock-ins. The new guarantee value is given by 

Gi+i = Gi+p{Gi+i-Gi)+ 

Between two lock in, the guaranteed amount is constant. 

We denote by Hjj_^ = G/+i/G/ the value of threshold immediately before the profit lock 

in. Similarly, we define XJ^^ = Ci^i/ HJ^^. In term of Xi and ^T+i guaranteed amount 
is rescaled by a factor 

fi+i{Xj,X-_^,) = ^ = l+p^Xy^^e -Xje '^'w' j 

2. Continuous lock-in. On rebalancing dates, some proportion p of the highest CPPI value attained 
is guaranteed. The new guarantee is (in fact, in this case, /(/) = / since lock-ins happen at 
every rebalancing dates) 

=max(G/, pG/+i) 

The multiplicative factor in this case is a function of XY_^_-^ alone and is given by 

fl^,{XT^,) = max (l, pXK.ie-^^.i^^^-^'^*) 

To get some intuition on what happens between ^^+1 ^^'-^ iii the case Xj > 1, see Figure [ij 

The scale is in unit of e . For periodic profit lock-in, there are three different regimes 

based on whether one started close to threshold or far above, the last regime being merely a transition 
between the first two. In the case of continuous profit lock- in, there is only one regime which is similar 
to the transition regime of the periodic case. 

4.2.2 Pricing 

The terminal payoffs we naturally want to consider are 

• Cn = GnXn, the CPPI strategy; 

• max(Gn,G„) = Gn max(Xn,l), the Guaranteed CPPI strategy; 

• {Gn — KGn)^ = Gn {Xn — K)'^ , a Call on CPPI with strike in percentage of the final guarantee; 
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3 4 

X before profit lock-in 



Figure 1: X/^i = X^^^/ fi{Xi^Xjj^^ as a function of Xj^^ in the case of a periodic profit lock-in 
with p = 0.5 for three values of u = Xje^ -^i r+sndt . ^ _ u = 2 and u = 2.85. The threshold 

value at I is equal to 0.95. 

• {KGn — Cn)^ = Gn {K — a Put on CPPI with strike in percentage of the final guarantee; 

• ^Cn<KG„i ^ digital Put with strike in percentage of the final guarantee {K = 1 gives the 
defeasance probability) ; 

• {Cn — K)^ , a Call on CPPI with strike in monetary unit; 

• {K — Cn)~^ , a Put on CPPI with strike in monetary unit; 

• ^c„<KG„i a digital Put with strike in monetary unit. 

With the exception of options with strike in monetary unit, these payoffs are of the form or 
P(Xn)- For a payoff which depends directly on the monetary value Ct the computation is different 



although the basic building blocks remain the same (see section 4.2.3). 

All computations throughout this section will be made with the hypothesis that there are no 
coupons. We will see at the end of the section that there is no difficulty in including them. 

With profit lock-in, the single variable Ci is no longer Markovian. One should introduce in addition 
the guarantee to construct a two-dimensional Markov model. However, the process ^ (where the 
threshold Hi is proportional to the guarantee Gi) is Markovian: the evolution of this quantity over one 
lock-in period depends only on its value at the beginning of the period, and periods are independent 
of each other. We introduce therefore two pricing kernels: 

/CW(x,y) = E[6{Xn-y)\Xo = x] 



/cW(x,y) = E 



—5{Xn - y) Xo 



(Without profit lock-in, both kernels are equal, this is the kernel /C of section 3.3 ) K,^^^ is used to 
compute the defeasance probability or digital options: 

E[P(X„)|Xo] = (/cWp)(Xo) (5) 



12 



whereas IC^^'^ is used to compute payoffs of the form GnP{X, 

Gn 



E[G„P(X„)|Xo] =GoE 



Go (/C(i)p) (Xo) 



(6) 



To compute the pricing kernels, the first step is to compute the evolution operator between two 
lock-in dates. We consider a process which ends just before the lock-in at date tz+i, that we denote 
by ^7+1- From section 3.3 the evolution operator on this period is given by 



^1,1+1- — ■ ■ ■ -^i(7+l)-l 



(7) 



Once this operator is computed, the lock-in has to be incorporated. It is given by a map which acts 
as 

x; 



X 



1 



fi+i{Xi,Xj^-^^) 

Compounding this map with the operator without profit lock-in gives the transition kernel on the 
lock-in period, including lock-in: 

vP {x,y)=E[ 5{Xj+i -y)\Xj = x]= [ dz ICjj+,- {x,z)6(y- 

When the multiplying factor of the guarantee is included, it gives the second kernel 

vP(x,y)=E ^^^^ 5{Xi+i - y)\Xi = x = /" dz/Cj j+i- (x, z) 5 ( y - - — j -]fi+i{x,z) 



As one can see in Figure [T| there might be, depending on the value of zero, one, two or an 

infinity of values of -Y^+i corresponding to different states of the portfolio before rebalancing. 

As the processes for Xj on different periods are independent, the total pricing kernels are given 
by the operator products 



They can be used in formulas ^ and ^ to compute the expected value of the quantity considered. 

To account for coupons, the only modification to make is that in the definition of the operator 
lCj jj_i- , Mi has to be replaced by Mi the evolution operator in the presence of coupons. The pricing 
then goes along the same lines as in the coupon case, making use of the same semi-direct product. 

4.2.3 Options v^ith strike in absolute amount 

With a strike in absolute amount, or equivalently in percentage of the initial guarantee, the compu- 
tation scheme is different. 

Here we write a transition matrix directly for the CPPI index value Cj. Technically we consider 
a deterministic rescaling which defines the variable 



ry Ci f^tr+SH)dt CiGi 

Zi = ~ 



Xi 



G, 



HiGo Gq 

We want to compute the transition operator for this variable 

Uix,y)=E[6{Zn-y) \ Zq = x] 
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This transition operator gives the probabihty distribution at the maturity of the CPPI and aUows 
the pricing of payoff P{Cn) as 

E[P{Cn)\Zo] = {UP){Zo) 

The transition operator can be computed by recursion. Introducing the evolution operator from 
lock-in / to maturity 



Ui{x,y)=E 
the basic step is the recursion formula 



Ui{x,y)= [ dzlCj j+i-{x,z)Ui+i( - — J -,- — J -) ] 

J \fi+i[x,z) fi+i{x,z)J fi+i{x,z) 

where }CJ J_^_l- is given by equation ([7]). This equation states that the evolution after a lock- in is the 
same evolution as without any lock-in but for a rescaled variable. The cumulative operator 

/y 
dzUi{x,z) 
-oo 

has a simpler recursion formula, more adapted to numerical interpolation: 

Qj{x,y)= / dz/C/ j+i-(x,z) Q/+i( - — ^ -,- — j -) 

J \fi+i{x,z) fi+i{x,z)J 

Using this equation backward from the CPPI maturity, the full transition matrix U = Uq can be 
computed and gives the distribution of the final CPPI index value, the price of the CPPI strategy or 
any option with absolute strike. Note that this computation scheme is no longer a Markov chain: it 
can not be reduced to a point wise rescaling followed by a matrix multiplication. 

4.3 Conditional rebalancing 

In order to minimize transaction costs and to simplify the management of the portfolio, one might 
choose not to rebalance the position if it is close enough to the previous one. For instance, one starts 
at time to at a given level Xq with the risky asset weighting level Wq = w{Xq). The CPPI will be 
rebalance on ti only if (the theoretical weighting) w{Xi) is such that w{Xi) ^ [(1 — f5)Wo, (1 -|-/3)Wo] 
for a given /? > 0. If w{Xi) lies within the interval, the risky asset weighting for the next period is 
left at Wq. This makes Wi a function of both Xq and Xi. On the following rebalancing date, W2 
would be a function of Xq, Xi and X2 and so on... To preserve the Markov property, we note that in 
this setup Wi is a function of Wj_i and Xi only: we have to follow the diffusion of both X and W , 
turning our problem 2 dimensional. The dynamics of W is given by 

Wi+i = Wil\ ^(x,+^) 

Instead of keeping track of VF, it is more convenient to change variable and set 

l+Y, 

The new variable Yi measures the the difference between the exposure Wi at time ti and its " theoretical 
value" w[Xi). If one defines 

w{Xi+i) ^ w{Xi+i) 
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the dynamics of Y is given by 



The main advantage of considering this new variable over direct propagation of W is that for 
typical rebalancing rules, /? is equal to a few percents and hence Y is constrained into a small 
interval which we can approximate with a small number of points in the Y-direction keeping the two 
dimensional problem tractable on a computer. The equation driving the dynamics of the CPPI over 
one period (see ([s]) and following) still hold if one substitutes w{Xi) with Wi. 

4.4 Open-ended CPPI 

An open-ended CPPI has no fixed maturity: the investor can exit at any time. It is therefore an 
American-style product. In addition, it has continuous profit lock-in: the guarantee is set to some 
proportion p of the highest value attained. As there is no fixed maturity, the threshold cannot be 
discounted, it is equal to the guarantee. 

As there the threshold is flat between rebalancing, the threshold spread is set over each period to 
the opposite of the interest rate sh = —r such that sh + r = 0. The continuous lock-in is handled as 
explained in section [4?2| on each period the evolution matrix Mi is used to compute the two pricing 
kernels vj''^^ and v| ^ The American feature is taken into account by propagating backwards the 
payoff vector, taking the maximum with the exercise payoff vector at each step. At date this gives 
the expected value of the American product conditional to the starting value of the CPPI index. 

If there is no terminal maturity, we force some long range maturity by hand, so that the process 
can be bootstrapped. If the put price diverges when this terminal maturity goes to infinity, the fair 
price of the product in the model considered is infinite, which probably means that the fees are not 
high enough to compensate for the gap risk. 



5 Numerical aspects 

In order to make practical computation, one discretizes the state space by a finite space (the grid) 
g = {gi, . . . , gN}- We also introduce an "separation" grid: g' = {g[, . . . , g'j^} such that gk-i < g'j^ < gk- 
Linear operators are approximated by matrices on this grid. Non linear features such as profit lock-ins 
are obtained by point wise operations on matrix element. Hereafter, we will denote the discretized 
operator corresponding to an operator O by O. In our implementation, the discretization scheme 

goes as follows: one replaces the probability density P [-'^j+i = y \ Xi = x\ by a matrix P ' whose 

elements -P^/ are given hy {1 < k,l < N) 

Pkt^'^ = I dyF = y I X, = <7fc] 

where The Markov property 

P [X,+2 = z\Xi = x]= J dyF [Xi+2 = z I Xi+i = y] F =y\Xi = x] 



is then replaced by the relation: 



N 

-p{i,i+2) ^ (i+l,i+2)p{i,i+l) 
^ km / J ^ Im ^ kl 



1=1 



This discretization scheme amounts to considering that, starting form a point Xi = x, if we arrive at 
a point Xi^i such that g'f^ < Xj+i < g'i^_^_i, the next step will start at X^+i = gi^. 
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5.1 Grids 



We take a fixed grid for the whole life of the product. The pricing could be improved by changing 
the grid extension as the maturity grows up, but it does not appear to be so important as long as 
the variance of individual steps is controlled. 

The upper bound of the grid is computed using the lognormal, continuous approximation. The 
process C — H is supposed to be log-normally distributed with volatility 'U'max'^ where tfmax is the 
maximum exposure, equal to the multiplier if no specific maximum exposure is set. a is the At-The- 
Money volatility of the underlying asset with the same maturity. Inverting the cumulative law of this 
lognormal distribution for some value 1 — e close to 1 provides an upper bound for the grid. Because 
of the leverage, this bound can be very high for a long maturity. 

The lower bound of the grid is constructed from the upper bound. If the risky asset goes to zero 
or close to zero, the CPPI index goes from C to (1 — w)C. Therefore the lower bound is chosen as 
(1 — Wmax) times the upper bound. It is generally a deeply negative value. 

Between these two points, the grid is densified around the threshold, where the behavior of the 
index changes and where the distribution can be sharp. If there is some cushion limit which stops 
investment in the risky asset (see section 7.5 for a precise definition), the grid is exponential between 
the threshold and the cushion limit: in this case there is only a drift motion in this region and a grid 
too large would introduce diffusion. We detail below how this diffusion is prevented using negative 
probabilities. 

The "separation" grid is introduced to separate points of the first grid and defines a partition of 
the real axis. When computing a transition probability, the probability to be at a given grid point 
is defined as the probability to lie inside the interval around this point, computed by difference from 
the cumulative function. This ensures that the sum of all probabilities is always 1. 



5.2 Matrix-Matrix vs Matrix- Vector multiplications 

When there is no explicit time dependence, the computation of the price amounts to propagating 
a vector through the n*^ of a matrix. Numerically, there are two ways of doing the computation. 
The first solution is to take the v}^ power of the matrix M by fast exponentiation |Alb07j : M^'' are 
computed recursively for growing p hy M^' = M^'-' M^'-\ The power n is decomposed in basis 2 
as 

i 

and M" is obtained by multiplying the corresponding M^''' : 

i 

For a discretization on a grid of size N ^ matrix multiplication has complexity 0{N^). The number 
of multiplication to perform is O(lnn). Thus this solution has complexity 0(ln(n)A'^^). The other 
solution is to sequentially apply n times the operator M to the vector P. The matrix vector multipli- 
cation being of complexity 0{N'^) this gives a total complexity in 0{nN'^). Depending of N , n and of 
the architecture used to perform the linear algebra, one will choose the fastest one for his particular 
needs. 



5.3 Mean and Variance control 

The discretization of the transition kernel introduce a bias in the mean value of the index C. To 
prevent this and get the correct forward price, the probability distribution constructed is slightly 
shifted on the grid in the opposite direction, such that the mean value computed on the grid equals 
the theoretical mean value, conditionally to the starting point. 
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Similarly, the variance of the process conditionally to the starting point is also adjusted. If the 
variance computed on the grid is too high, the probabilities above and below the mean value are 
slightly displaced closer to the mean value, in a way which preserves the mean. The converse is done 
when the variance computed on the grid is lower than the theoretical variance. 

Finally, we get probability distribution conditional to the starting point which gives the correct 
mean and variance computed on the grid. This appears to be crucial for convergence when the 
number of rebalancing dates grows up. If the mean is not controlled, the expected value of the CPPI 
strategy can diverge. If the variance is not controlled at each step, the total variance is not correct 
and a reasonable estimate of the gap probability can require a large number of points on the grid, 
which slows down the computation in N"^. 



5.4 Convective limit 

When the threshold does not grow at the same rate as the risk-free rate, the process X has some 
drift. In region where the pure diffusive part is small (i.e. near the threshold), this drift might 
have undesired effects: unless the drift over one rebalancing period exactly corresponds to an integer 
number of steps in the grid, the discretization scheme will introduce some spurious diffusion. This 
enlarges the distribution in this important part of the grid where the threshold lives and can lead 
to very inaccurate results for the gap probability and mean value. This phenomenon is more easily 
identified when the CPPI includes a cushion limit in which case the CPPI structure is fully monetized 
before reaching the threshold. 

An efficient way to suppress this problem is to allow negative probabilities such that the variance 
can be set to zero even when the process must be between two points of the grid. Let consider that 
starting from the grid point gi, the process must drift to gi-i < x < gi. If no negative probabilities 
are allowed, the solution which minimizes the variance while preserving the mean value is to take 
probabilities 

X - Qi 



Pi-i = 

9i-i - 9i 

Pi = 1- pi_i 

The mean value is x but the variance is not zero. We consider instead the following "transition 
probabilities" : 

Pi-3 = aq 

Pi-2 = -q 

Pi~i = P 

Pi = l-p+{2-a)q 

Pi+i = -q 

where p and q are positive and a is a numerical constant. Taking a non vanishing a helps stabilizing 
the process. We set it to a = 0.1 as this value gives good results for both stability of the process and 

convergence of the gap estimators, p and q are chosen so that the mean value is x and the variance 
vanishes. These conditions translates into a linear system of two equations in p and q which is easily 
solved. We bound the negative probability q to 0.05 to avoid explosions of negative probabilities in 
some parts of the grid. 

This drift scheme appears to be impressively better than the naive one. The accuracy obtained 
for a grid of 2000 points can be reached with a few hundreds of points. Drawbacks of negative gap 
probabilities or gap mean occur only for very sparse grids. Further work would be necessary to be 
able to control a priori the stability of this trick. 
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5.5 Semi-direct product in matrix form 



The semi-direct product we defined in the section 4.3 is easily handled with matrices. If N is the size 
of the discretization, to a pair (O, /) where O is an operator and / is a function, we associate the 
discretization given by the (N + 1) x {N + 1) block matrix 



V" 



o f 



\ 



1 / 

The product (Oi, /i) {O2, f2) of two such matrices is given by 

/ 





> 


Oi O2 







1 



which is precisely the same algebra as the semi-direct product we previously defined. 

The action of a pair (O, /) on a function g is now represented by the standard product of matrices 

/ 

V 





\ 


/ \ 


/ \ 





/ 


g 


O-g + f 





1 / 


\ 1 / 


\ 1 / 



5.6 Piecewise constant approximation 

Quantities involved in the computation of a single period transition matrix depend on the precise 
features of the CPPI. There is at least the variance of the risky asset over the considered period. When 
spreads over the risk-free rate are applied, the value of the zero-coupon spread over the period can 
change, if the length of the period changes a bit, depending on the day count basis, or if the spread 
is not expressed in continuous mode and the interest rate is not constant. Moreover when there 
are coupons or continuous profit lock-in, the zero-coupon value between the rebalancing date and 
the maturity of the CPPI enters also the computation. In these cases, an exact computation would 
require computing all single period matrices. The computation can be fastened if such quantities are 
approximated by piecewise constant functions. Technically, we choose a tolerance for the variation 
of variables and we compute only one matrices for periods with parameters which do not differ by 
values larger than this tolerance. Mean values for the parameters are computing in the following way. 
The variance of the risky asset is the arithmetic mean of the variance over single periods whereas 
for zero coupons it is the geometric mean which is used, such that their product is equal to the 
real zero-coupon over the sum of periods considered. The numerical impact of this approximation is 
studied in section ID 



6 Numerical convergence 

In order to analyze the impact of the grid size, we set up two benchmark CPPIs and study the price 
of a put option with a strike equal to the guarantee. Both CPPI start on November 12, 2008 and 
end on November 7, 2018, their multiplier is 4, their nominal is 1 000 000 and are rebalanced every 
week. On the start date, the risky underlying is quoted at 3207. 

• Vanilla CPPI. The rate curve is flat with a constant rate of 5%, the volatility of the risky 
underlying is 50%. Threshold is natural and there are no fees. 
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• Non- Vanilla CPPI. Volatility is still 50% but the rates are time dependent. The threshold 
is linear with an initial value of 60% of the nominal. There are a cushion limilQ of 5% and 
proportional fees of 50bp. 

In the following, we consider a European put option with a strike given by the guaranteed nominal 
(there is no profit lock- in). The computation date is chosen to be November 16, 2008 and the spot 
of the risky underlying on that date is 3190. 

The computations (both Markov and Monte Carlo methods) were performed on a PC equipped 
with an Intel Core2 6700 at 2.66 GHz and 3.25 GB of memory. For the Markov method, we used 
Intel's Math Kernel Library to perform matrix-matrix and matrix-vector multiplications. The source 
codes are not heavily optimized and hence computing times are only relevant for comparison purposes. 
In the case of Markov, the times we indicate correspond to only one pricing which is sufficient for 
getting the price, the delta and the gamma of the option. In order to get the vega two additional 
pricing are performed (but only one more could be enough). 



6.1 Vanilla CPPI 

In the case of the vanilla CPPI, we can compute analytically the price and the sensibilities of the 
option using the formulas of appendix|B| The results obtained are (these numbers have to be compared 
to the CPPI initial value of 1 000 000): 



Put Price 170.5530 

Put Delta 0.2177 

Put Vega 68.2553 

Gap Proportion (%) 9.7989 

In table [T] we give the relative errors with respect to the analytical results and the computation 
time, depending on the grid size. Even with a small number of points, there are almost no pricing 
errors: this is due to the fact that we (try to) force the distributions on each rebalancing period to 
have the correct mean, variance, conditional (on defeasance on that rebalancing date) mean and con- 
ditional variance: the price of the option and the gap proportion of the CPPI are then automatically 
exact. The vega being computed by finite difference (with a fixed bump) on the price has almost no 



Grid Size 


Price 


Delt 


a 




Vega 


Gap Prop. 


Time (s) 


50 


-3.09 X 10"*^ 


-6.57 X 


10- 


-3 


-6.34 X 


10" 


-6 


-3.32 X 10" 


10 


0.06 


100 


-9.02 X 10-^ 


-6.89 X 


10- 


-5 


-6.52 X 


10- 


-6 


-3.32 X 10- 


10 


0.07 


250 


-4.50 X 10"^ 


-2.68 X 


10" 


-6 


-6.20 X 


10- 


-6 


-3.32 X 10- 


10 


0.10 


500 


-5.58 X 10-8 


-6.35 X 


10- 


-7 


-6.42 X 


10" 


-6 


-3.32 X 10- 


10 


0.24 


750 


-5.38 X 10-8 


-3.57 X 


10" 


-7 


-6.41 X 


10" 


-6 


-3.33 X 10- 


10 


0.59 


1000 


-5.68 X 10-8 


-2.55 X 


10- 


-7 


-6.42 X 


10" 


~6 


-3.33 X 10- 


10 


1.49 


1500 


-5.71 X 10-8 


—1.66 X 


10- 


-7 


-6.42 X 


10" 


-6 


-3.33 X 10- 


10 


3.80 


2000 


-5.67 X 10-8 


-1.31 X 


10" 


-7 


-6.42 X 


10" 


-6 


-3.32 X 10- 


10 


7.33 


2500 


-5.69 X 10-8 


-1.09 X 


10" 


-7 


-6.43 X 


10" 


-6 


-3.32 X 10- 


10 


12.41 



Table 1: Relative errors for the put on the vanilla CPPI. 

dependance on the grid size. The delta is computed by finite difference on the grid and hence slightly 
depends on the grid size by a third order effect. The non- vanilla case will give more realistic errors 
and provide a better estimate of the necessary grid size and computation time. This test case was 
intended to check the algorithm in a case where closed formulas were known. 



See section 7.5 
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6.2 Monte-Carlo pricing 



We compare our results to a standard Monte-Carlo pricing run on the same computer. We take a 
Sobol random generator and transform the uniform laws into normal laws through a Moro algorithm. 
In the continuous rebalancing limit, the CPPI terminal value depends only on the underlying spot at 
maturity. So we use a Brownian bridge path generation which takes advantage of this property and 
improves convergence. Figure [2] shows the convergence graph up to 20 000 000 paths. In 18 hours, 
the price and greeks have not convergecPj 



6.3 Grid size 

To study the influence of the grid size on the outputs we consider now the non-vanilla case. Table 
|2] shows the put price, delta, gamma, and defeasance probability depending on the grid size. The 
precision is correct from 100 points, and quite good from 500 points. Note that the goal of this 
example is to study the numerical precision of the pricing model: as the vega is very high compared 
to the price, even with a grid size of 100 the numerical errors are negligible in front of model risk. 



Grid Size 


Price 


Delta 


Vega 


Gap Prop. 


50 


2028.6075 


3.3718 


-48.0019 


-6.7264 X 10" 


-3 


100 


129.7479 


0.1902 


52.9596 


5.2489 X 10~ 


-4 


250 


133.4823 


0.1941 


54.4288 


5.8667 X 10" 


-4 


500 


138.3577 


0.2016 


56.3636 


5.9074 X 10" 


-4 


750 


139.3600 


0.2034 


56.7991 


5.7116 X 10" 


-4 


1000 


139.7006 


0.2040 


56.9460 


5.6653 X 10" 


-4 


1500 


139.9571 


0.2044 


57.0554 


5.6465 X 10" 


-4 


2000 


140.0305 


0.2045 


57.0868 


5.6360 X 10^ 


-4 


2500 


140.0549 


0.2046 


57.0972 


5.6350 X 10- 


-4 



Table 2: Convergence results for the put on the non-vanilla CPPI. 



6.4 Piecewise constant approximation 

In order to speed up the pricing, the interest and threshold rates are approximated by piecewise 
constant functions as well as variance of the underlying between two rebalancing dates. The approxi- 
mation is controlled by a tolerance number: when the rate has moved by more than this number a new 
constant is used. Constant rates are calibrated so that zero-coupons are exact on the period. This 
reduces in a large proportion the number of matrices to be generated. Computation times depending 
on the tolerance level and the grid size are presented in table |3j Relative errors with respect to the 
most accurate case (exact rates, 2500 points) are presented in tables |4] [5| |6] and [T] A good precision 
can be reached in about Is. 



^ Other variance reduction techniques could have been used. With importance sampling where more weight would 
be given to tail events. The knowledge of the closed formula in the vanilla case can also be used to implement a 
control-variate for non-vanilla cases. However jumps in the price due to extreme events would not necessarily occur at 
the same time for both products: this could produce a very jumpy estimate of the covariance. It is not clear if the 
convergence would be better. 
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Figure 2: Convergence of the Monte Carlo method for price, delta and vega in the vanilla case, as 
functions of the number of paths. In each graph, the dashed line indicates the corresponding exact 
value obtained with analytical computation. 
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-3 
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-3 
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2 


4.00x10-2 


1500 


109.31 




10.20 




5.69 


4.41 




4.17 




-6.99x10- 


-4 


-5.11x10- 


-4 


9.28x10--* 


1.24x10- 


2 


4.15x10-2 


2000 


173.11 




17.70 




10.46 


8.46 




8.17 




-1.75x10- 


-4 


1.47x10- 




1.46x10-3 


1.29x10- 
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4.22x10-2 


2500 


266.14 




28.97 




18.06 


14.99 




14.46 








1.89x10- 


4 


1.62x10-3 


1.31x10- 


2 


4.24x10-2 



Table 3: Computing time in seconds as a function of the approximation level and the grid size. The 
numbers below indicate the relative error on the non-vanilla put price with respect to the best estimate. 
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3.83 X 10" 


-2 


1000 


-2.53 X 10-3 


-2.34 X 


10-3 


-8.90 X 10-4 
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4.00 X 10" 
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-6.99 X 10-4 


-5.11 X 
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9.28 X 10-4 


1.24 X 10-2 


4.15 X 10" 


-2 


2000 


-1.75 X 10-4 


1.47 X 
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1.46 X 10-3 


1.29 X 10-2 


4.22 X 10" 


-2 
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1.89 X 
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1.62 X 10-3 


1.31 X 10-2 


4.24 X 10" 


-2 



Table 4: Relative errors for the price of the put on the non-vanilla CPPI as a function of the approx- 
imation level and the grid size. 





0% 


5% 






20% 


50% 




100% 




50 


1.55 X 10+1 


1.55 X 


10+1 


1.56 X 10+1 


1.56 X 


10+1 


1.62 X 10+1 


100 


-6.99 X 10-2 


-6.98 X 


10- 


-2 


-6.91 X 10-2 


-6.26 X 


10- 


-2 


-5.27 X 10- 


-2 


250 


-5.10 X 10-2 


-5.09 X 


10- 


-2 


-5.04 X 10-2 


-4.51 X 


10- 


-2 


-3.47 X 10- 


-2 


500 


-1.44 X 10-2 


-1.43 X 


10" 


-2 


-1.38 X 10-2 


-8.53 X 


10" 


-3 


1.30 X 10- 


-3 


750 


-5.54 X 10-3 


-5.48 X 


10- 


-3 


-4.95 X 10-3 


1.40 X 


10- 


-4 


1.03 X 10- 


-2 


1000 


-2.76 X 10-3 


-2.69 X 


10- 


-3 


-2.16 X 10-3 


3.01 X 


10- 


-3 


1.34 X 10" 


-2 


1500 


-9.01 X 10-4 


-8.33 X 


10- 


-4 


-2.94 X 10-4 


4.92 X 


10" 


-3 


1.54 X 10- 


-2 


2000 


-2.85 X 10-4 


-2.17 X 


10- 


-4 


3.21 X 10-4 


5.53 X 


10- 


-3 


1.60 X 10- 


-2 


2500 




6.82 X 


10- 


-5 


6.12 X 10-4 


5.83 X 


10- 


-3 


1.63 X 10- 


-2 



Table 5: Relative errors for the delta of the put on the non-vanilla CPPI as a function of the approx- 
imation level and the grid size. 
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0% 




5% 




20% 




50% 




100% 




50 


-1.84 X 10+" 


-1.84 X 10"^" 


-1.79 X 10+" 


1 Q1 V 




hU 


-2.32 X 10+" 


100 


—7 25 X 10" 


-2 


-7 23 X 10" 


-2 


-7.11 X 10" 


-2 


1 9 V 


1 0" 


-2 


-3 35 X 10" 


-2 


250 


—4.67 X 10" 


-2 


—4 66 X 1 0" 


-2 


—4 53 X 10" 


-2 


— O.Oa A 


1 n~ 


-2 


—1 00 X 10" 


-2 


500 


-1.28 X 10" 


-2 


-1.27 X 10" 


-2 


-1.14 X 10" 


-2 


— 1 . / tj A 


1 0" 


-3 


2.58 X 10" 


-2 


750 


-5.22 X 10" 


-3 


-5.05 X 10" 


-3 


-3.76 X 10" 


-3 


6.33 X 


10" 


-3 


3.35 X 10" 


-2 


1000 


-2.65 X 10" 


-3 


-2.48 X 10" 


-3 


-1.20 X 10" 


-3 


8.74 X 


10" 


-3 


3.55 X 10" 


-2 


1500 


-7.32 X 10" 


-4 


-5.66 X 10" 


-4 


7.03 X 10" 


-4 


1.06 X 


10" 


-2 


3.72 X 10" 


-2 


2000 


-1.82 X 10" 


-4 


-1.44 X 10" 


-5 


1.26 X 10" 


-3 


1.11 X 


10" 


-2 


3.79 X 10" 


-2 


2500 






1.66 X 10" 


-4 


1.43 X 10" 


-3 


1.13 X 


10" 


-2 


3.80 X 10" 


~2 



Table 6: Relative errors for the vega of the put on the non-vanilla CPPI as a function of the approx- 
imation level and the grid size. 





0% 




5% 




20% 




50% 




100% 


50 


-1.29 X 10^ 


-1 


-1.29 X 10^ 


-1 


-1.29 X 10^ 


-1 


-1.25 X 10+1 


-9.52 X 10+° 


100 


-6.85 X 10" 


-2 


-6.48 X 10" 


-2 


-4.92 X 10" 


-2 


7.58 X 10" 


-3 


7.34 X 10-2 


250 


4.11 X 10" 


-2 


4.27 X 10" 


-2 


5.47 X 10" 


-2 


1.14 X 10" 


-1 


3.44 X 10-1 


500 


4.83 X 10" 


-2 


5.05 X 10" 


-2 


6.35 X 10" 


-2 


1.28 X 10" 


-1 


4.43 X 10^1 


750 


1.36 X 10" 


-2 


1.50 X 10" 


-2 


2.56 X 10" 


-2 


1.07 X 10" 


-1 


3.72 X 10-1 


1000 


5.37 X 10" 


-3 


6.02 X 10" 


-3 


1.67 X 10" 


-2 


8.96 X 10" 


-2 


3.21 X 10^1 


1500 


2.04 X 10" 


-3 


3.67 X 10" 


-3 


1.28 X 10" 


-2 


8.19 X 10" 


-2 


3.01 X 10-1 


2000 


1.67 X 10" 


-4 


1.12 X 10" 


-3 


1.24 X 10" 


'2 


8.17 X 10" 


-2 


3.08 X 10-1 


2500 






2.41 X 10" 


-3 


9.33 X 10" 


-3 


7.84 X 10" 


-2 


3.06 X 10-2 



Table 7: Relative errors for the gap proportion of the non-vanilla CPPI as a function of the approxi- 
mation level and the grid size. 

For comparison, the convergence graphs of the Monte-Carlo pricing are shown on figure |3j As 
in the vanilla case, the price and greeks have not converged after 20 000 000 paths and 18 hours of 
computation. 

6.5 Profit lock-in 

Computing times are higher when profit lock-in is used. In this case, matrix-matrix multiplications 
have to be performed instead of matrix-vector multiplications. In addition, transition matrices over 
lock-in period have to be transformed, which can take some time. Both aspects make the algorithm 
0{N^) instead of 0{N'^) where is the grid size. For example, the N = 500 case with 20% tolerance 
on the rate variation takes 26.7 s instead of 510 ms. However this is still reasonable and much better 
than a Monte-Carlo pricing. Parallelization of the algorithm would also reduce the computation time. 



23 




100 1000 10000 100000 lc+06 lc+07 

Sampling Count 




100 1000 10000 100000 lc+06 lc+07 

Sampling Count 




_20 I • • ^ • • ^1 • • ^ • • ^ • • ^ 

100 1000 10000 100000 lc+06 lc+07 

Sampling Count 



Figure 3: Convergence of the Monte Carlo method for price, delta and vega in the non-vanilla case, 
as functions of the number of paths. In each graph, the dashed line indicates the most accurate value 
obtained with the Markov scheme. 
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7 Features analysis 



In order to analyze various features of CPPI strategies, we consider now a 10 years CPPI with 
weekly rebalancing on some risky asset with 20% volatility. The multiplier is 4 and there are no 
fees nor spreads on the risk-free asset. The risk-free interest is taken constant at 5% in continuous 
compounding. In a perspective of pricing and hedging the gap risk, we consider the risk-neutral 
probability. Other probability spaces could be used for other purposes. 
For each case, we compute the following gap indicators 

• Gap Proportion: probability of the CPPI strategy being below the guarantee at maturity, 
P[Ct < 1]; 

• Conditional Loss: mean value of the loss when the strategy ends below the guarantee, E[(l — 
Ct)^\Ct < 1]; 

• Expected Loss: expected value of the loss without conditioning, E[(l — Ct)^]- 
7.1 Risky asset model 

The model of the risky asset has a large influence on the gap risk. We have performed computations 
for three cases: 

• Brownian motion with 20% constant volatility; 

• Jump-diffusion, adding exponentially decreasing jumps to the Brownian motion in both direc- 
tions (Kou model). The volatility of the Brownian is still a = 20% and the jumps param- 
eter set is A+ = 0.1, ?7+ = 0.05, A_ = 0.1, r/_, = 0.1. The total volatility of the process is 
y^cT^ + 2X+ril + 2A_r/2 ~ 21.6795%. 

• Brownian motion which gives the same expected loss E[(l — Ct)~^] as the jump diffusion case. 
The volatility which fits this number is close to 62%. 

As shown in table [s] and figure |4| adding jumps changes completely the loss. With a 20% volatility 
in pure Brownian motion, there is no loss probability up to numerical accuracy (~ 10"^'^). With 
multiplier 4, a loss occurs when the risky asset price drops by more than 25%, which have negligible 
probability in one week at the volatility considered. Adding jumps modifies completely this behavior 
and results in a completely different situation, more consistent with observed skews. Extreme moves 
are allowed with some small probability. With the parameters that we have taken, this gives a 5% 
probability of breaking the threshold in ten years, which is not negligible. 



Model 


Gap Prop. 


Cond. Loss 


Expect. Loss 


Brownian 20 % 
Kou 20 % 
Brownian 62 % 


0.00 % 
5.71 % 
21.78% 


0.00 % 
18.41 % 
4.85 % 


0.000 % 
1.052% 
1.058% 



Table 8: Gap proportion, conditional loss and expected loss depending on the risky asset model. 

In an attempt to take into account the gap risk using pure Brownian motion, one could be tempted 
to increase the model volatility. We find that the volatility which would give the same expected loss 
would be as high as 62%. However it results in a completely different probability distribution for the 
strategy price. The loss are much more probable but have a much lower mean value: the left tails 
have different shapes. The right part of the graph depends mainly on the variance of the underlying 
process and is almost given by a lognormal law, shifted by 1, centered on 1.65 (the forward value with 
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Figure 4: Probability Density of a CPPI strategy final value depending on the risky asset model. CPPI 
has maturity 10 years, weekly rebalancing, multiplier 4- The initial investment and the guaranteed 
level are 1. The risk- free interest rate is 5%, so the mean price of the strategy is 1.65. There are 
no fees. The first curve is the probability density for a risky asset following a Brownian motion 
with 20% constant volatility, for the second curves, jumps have been added: downward jump with 
intensity 0.1 and mean value 0.1, upward jumps with intensity 0.1 and mean value 0.05. The third 
curve corresponds to a Brownian motion without jumps which gives the same gap risk, i.e. with 62% 
volatility. 

5% interest rate) and with standard deviation mo\fT where m is the multipHer and o\/T the square 
root of the variance of the underlying. Therefore the probabihty density with 62% volatility has 
nothing to do with the 20% case neither in the gap part of the distribution nor in the well-performing 
part. Increasing the model volatility is not a solution for pricing or hedging the gap risk. 

For equities, jump- diffusion is more realistic model than pure Brownian motion: it is more consis- 
tent a non-vanishing short-term skew and historical returns seems to be consistent with Kou jump- 
diffusion model |CT07j . Therefore we consider in the following only this case, with the volatility and 
jump parameters given above. 

7.2 Multiplier 

The multiplier increases the variance of the CPPI strategy by leveraging the risky asset. As a 
consequence a big multiplier increases the gap risk and the probability of very high returns, as it can 
be seen on table |9] and figure |5j The case of unit multiplier is degenerate: there is no leverage and 
no need to rebalance. The strategy reduces to a zero-coupon plus some shares which are kept until 
the maturity. It can be checked that in this case there is no gap risk: the probability density is null 
below 1. 

The floor is broken when there is a downward move larger than the inverse of the multiplier 
between two rebalancing dates: ^ < — — . For a multiplier larger than 1, there is a positive 
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Figure 5: Probability Density of the CPPI strategy depending on the multiplier. The underlying is the 
Kou model which includes jumps with the same parameters as in figure The first graph shows the 
shape of curves near the threshold. The second graph shows the tails of the distributions. 
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IVTnltiDlipT" 


Oar) PvoD 


Oonrl TjOss 




1 


0.00% 


0.00 % 


0.000 % 


2 


0.10% 


5.91 % 


0.605 % 


4 


5.71 % 


18.41% 


1.052% 


6 


15.51% 


33.37% 


5.177% 



Table 9: Gap proportion, conditional loss and expected loss depending on the multiplier. 

probability of losing more than the initial investment. Because of the leverage, the CPPI strategy 
can end with a negative value when a downward jump in the risky asset is larger than the financing 
component. 

7.3 Maximum exposure 

In order to reduce the exposition to adverse moves and reduce the gap risk, an efficient way is setting 
a maximum exposure, which is a cap on the risky asset weighting. When there is no such constraint, 
the asymptotic value of the exposure is given by the multiplier. For the same CPPI as above, we 
show in table [To] the effect on the gap risk of reducing the maximum exposure from the unconstrained 
400% case (recall that m = 4 in our example). 



Max. Expo. 


Gap Prop. 


Cond. Loss 


Expect. Loss 


400% 


5.71 % 


18.41% 


1.052% 


200% 


5.16% 


8.13% 


0.419% 


150% 


4.44% 


5.88% 


0.261 % 


100% 


2.92% 


3.66% 


0.107% 



Table 10: Gap proportion, conditional loss and expected loss depending on the maximum exposure. 

When a loss occurs on the CPPI strategy, its value is proportional to the index level and the 
multiplier. Therefore bounding the exposure decreases the probability of a large loss by capping the 
exposure when the index level is high. The probability of little losses are less modified: the left tail 
is shortened. The effects on both tails of the distribution can be seen in figure [6] There is a slope 
change at the level where the maximum exposure is reached. The power effect of the CPPI strategy 
is stopped and the right part of the distribution is more or less a dilated version of the underlying 
probability density. The second graph shows that both tails are reduced; in particular limiting the 
exposure to 100% obviously prevent the CPPI strategy to drop below 0. 
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Figure 6: Probability Density of the CPPI strategy depending on the maximum exposure. The multi- 
plier is 4- The first graph shows the shape of curves near the threshold. The second graph shows the 
tails of the distributions. 
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7.4 Minimum exposure 



One may want to set a minimum exposure, i.e. a floor on the risky asset weighting, so that there 
is always some possibihty to recover from a gap event. Figure [7] shows that this features does not 
change the CPPI distribution above some value which is not very high, here 1.2, to be compared with 
the mean value which is 1.65. The gap risk increases when the minimum exposure gets higher, as 




Figure 7: Probability Density of the CPPI strategy depending on the minimum exposure. The tails 
are not shown as they are almost identical. 



shown in table 11 Some small diffusion losses appear, which results in a large gap proportion but 
a relatively low increase of the expected loss. Large losses come from jumps starting at high levels, 
where the distribution is not changed. 



Min. Expo. 


Gap Prop. 


Cond. Loss 


Expect. Loss 


0% 


5.71 % 


18.41 % 


1.052% 


5% 


19.77% 


6.39 % 


1.263% 


10% 


24.87% 


6.70 % 


1.666% 


20% 


30.80 % 


8.96 % 


2.759% 



Table 11: Gap proportion, conditional loss and expected loss depending on the minimum exposure. 



7.5 Cushion limit 

Capping the exposure to some maximum value reduces the gap risk. An other feature which can 
prevent some losses is the cushion limit: the exposure is set to zero below some cushion value. Gap 
indicators of table 12 show that this effectively reduces the gap proportion: the cases where the CPPI 
strategy fall below the limit value are frozen and cannot defease any longer. However the expected 
loss is not significantly reduced. 
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Limit 


Oar) PvoD 


Oond TjOss 


J l_/ V_> V_> LJ ft 1 /i 11 1 


0% 


5.71 % 


18.41 % 


1.052 % 


3% 


3.95 % 


26.22% 


1.036% 


8% 


2.92% 


33.83 % 


0.988 % 


15% 


2.05 % 


43.05 % 


0.884 % 



Table 12: Gap proportion, conditional loss and expected loss depending on the cushion limit. 



The effect of this feature on the distribution can be observed on figure [8j the Dirac-like probability 
peak just above the threshold is displaced because the motion is stopped before low values are reached. 
Losses due to large jumps from higher values are not modified. 



7.6 Profit lock-in 

The profit lock-in feature locks some performance at periodic intervals. We consider here annual 



lock-in and look at the effect of the lock-in proportion. Table 13 presents two series of indicators. In 
the first set of column, the guarantee is kept at the initial level whereas the threshold value used to 
compute the cushion is increased. The gap risk is therefore reduced because the lock-in only reduces 
the exposure. The second set of column deals with the case where the guaranteed level is increased 
of the lock-in amount. The gap proportion is not changed in this case, as the CPPI threshold and 
guarantee are increased in the same proportion. However the expected loss is reduced because the 
exposure is reduced. 





Initial Guarantee 


Locked-In Guarantee 


Prop. 


Gap Prop. 


Cond. Loss 


Exp. Loss 


Gap Prop. 


Cond. Loss 


Exp. Loss 


0% 


5.71 % 


18.41 % 


1.052% 


5.71 % 


18.41 % 


1.052% 


25% 


2.20% 


18.95% 


0.417% 


5.71 % 


13.61 % 


0.777% 


50% 


1.72% 


14.57% 


0.251 % 


5.71 % 


10.21% 


0.583 % 


75% 


1.52% 


12.91% 


0.197% 


5.71 % 


7.84 % 


0.448 % 



Table 13: Gap proportion, conditional loss and expected loss depending on the Lock-In proportion, 
with respect to the initial guarantee and with respect to the final locked-in guarantee. 

Figure |9] shows the effect on the curve and on its tails which are reduced by lowering the leverage. 
The effect is smoother than with a maximum exposure, but if the lock-in applied to the guarantee it 
does not reduce the gap risk as efficiently. 



7.7 Artificial threshold 

Finally we discuss the effect of artificial thresholds: the threshold is no longer defined by the natural 
zero-coupon value. Either it is defined by a linear function of time, by a fixed rate, or by a spread 
over the zero-coupon curve. We focus here on the latter case. With a negative spread, the initial 
threshold is higher than the natural one and finishes at the the level. The exposure is thus reduced, 
especially for low values of the CPPI index. This moves slightly a part of the distribution which was 
close to the threshold and smoothes the peak of the distribution. On the other hand, the right tail is 
also reduced as the mean value remain the same, here 1.65. 
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Figure 8: Probability Density of the CPPI strategy depending on the cushion limit. The second graph 
shows the left tail of the distributions. The curves are almost identical on the right tail. 
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Figure 9: Probability Density of the CPPI strategy depending on the profit lock-in proportion. Annual 
performance is locked with the given proportion. The second graph shows the tails of the distributions. 
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For positive spreads, the effect is tlie opposite: tlie tliresliold starts below its natural value and 
goes up with time to reach the same level at the CPPI maturity. As the CPPI distribution tends 
to concentrate just above the threshold value, this results in the threshold progressively encroaching 



some part of the distribution and produces large amounts of small losses as shown in figure 10 



This effect could be compensated by introducing a cushion limit. Thus the effect of increasing the 
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Figure 10: Probability Density of the CPPI strategy depending on the threshold spread. 

exposure at the beginning of the CPPI life would remain without this mechanical increase of the 
gap proportion. However increasing the exposure automatically increases the risk of large losses: the 
conditional loss increases also as shown in table [T4l 



Spread 


Gap Prop. 


Cond. Loss 


Expect. Loss 


-2% 


1.66% 


35.61 % 


0.592 % 


-1% 


2.47% 


31.52% 


0.779% 


0% 


5.71 % 


18.41 % 


1.052% 


0.5% 


48.16% 


3.87% 


1.862% 


1% 


59.33 % 


5.41 % 


3.212% 



Table 14: Gap proportion, conditional loss and expected loss depending on the threshold spread. 

When a linear threshold is chosen, depending on the initial level both regimes can exist succes- 
sively: the linear threshold can grow faster than the natural threshold and slower at the end. The 
effect of both regimes will partially compensate each other. 
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8 Conclusion 



Wc have introduced an efficient scheme for the pricing of CPPI strategics and options. Instead 
of foUowing the underlying spot process as in a Monte-Carlo pricing, only the CPPI strategy is 
considered. CPPIs are very sensitive to tails events; these are correctly handled numerically so that 
loss estimates can be computed for various underlying process. CPPI value distribution depends 
mainly on the mean volatility and on the left tail of the short term distribution for the underlying. 
With pure Brownian motion, huge volatility levels arc necessary to get some gap risk: deeply Out- 
The-Money forward volatility controls this gap risk. On the other side, the average evolution of 
a CPPI depends on average volatility until maturity. These two parameters should be controlled 
separately. Introducing jumps is the simplest way to handle this problem in a homogeneous model. 
We performed some analysis for exponentially distributed jumps, other jump distribution can be 
used as well. As the pricing depends on forward volatility, a direct and natural extension would 
be to introduce stochastic volatility through a few volatility regimes (with jumps depending on the 
volatility level) . This would be useful for CPPI management as the potential gap size increases with 
the leverage. 

The numerical analysis presented in this paper give some indications on the effect of various 
features aimed at reducing the gap risk. The most efficient features are the ones which reduce the 
leverage of the structure. Other techniques, such as setting a cushion limit (where the CPPI is 
completely de-leveraged before reaching the floor), only reduce the diffusive risk which is already 
almost null. 

We have not addressed the important question of hedging, which is beyond the scope of this 
work. The hedging of the regular part of the CPPI, where the threshold is not broken, is quite 
simple: the delta of the strategy is more or less equal to the risky asset weighting dictated by the 
strategy definition. The difficult part is gap risk hedging which requires instruments sensitive to large 
downward jumps of the underlying process. If such instruments exist, the pricing scheme described 
here can be used to compute hedging ratios for optimal hedging and get estimators of the residual 
risk. 

Acknowledgements 

We warmly thank Martial Millet for numerous stimulating discussions and for his constant availability. 
We also thank Erwan Curien for giving us the opportunity to work on this subject. Both of them are 
also acknowledged for a careful reading of this paper and useful comments. This work was supported 
by Sophis Technology. 



35 



Appendix A Diffusion models 

Black-Scholes 

If there are no dividends, under the T-forward risk neutral measure, the forward price of S at maturity 
is a martingale 

F 

This leads to the finite variations 



Fi 



With constant volatilities the cumulative and the (conditional) expected value of the variations 
of F up to some strike K are given by the usual formulas: 



Fi 



< K 



E 



Fi+i^ 



Af{d-) 



with 



d± 



In K±^a\ti+i-ti) 



<7^U+i — t 

where M denotes the cumulative function of a standard normal law 



Jump-diffusion 

Jumps can be added to the diffusion process in order to take into account large deviations in short 
timescale. The jump-diffusion model introduced by Kou [ Kou02j has three independent components: 
a Black-Scholes diffusion, downward jumps and upward jumps both exponentially distributed with 
distinct parameters. The asset price is the exponential of the sum of the three components. As a 
consequence, the distribution of the asset price logarithm is the convolution of three distributions. It 
can be efficiently computed by Fourier transforms. 

The diffusion is given by the following equation (again without dividends and under the T-forward 
risk neutral measure) 

dlnF = jdt + adW + dJ+ + d J_ 

dW is a standard Brownian motion, dJ± are jump processes with intensity X± and mean jump size 
r]±. The jumps' size distributions are 

^^±(^) = — e^^/^''±l±.>o 
V± 

The Brownian process dW, jump arrival times and jump sizes are all independent one from another. 
To ensure that T is a martingale, the parameter 7 must be set to 

7 = ^^A+ H A_ 

Moreover, the mean size of positive jumps 77+ must be lower than 1, i.e. the asset price cannot 
double or more during an infinitesimal period of time, in order to define a well behaved process]^ 

^Similarly, the variance of the process S is infinite for 77+ > | and, more generally, the n^^ moment will diverge if 
77+ > 
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Assuming constant coefficients for simplicity, tlie distribution 
Fourier transform from the characteristic exponent ip{u): 



of 



can be computed by 



with 



1 

2^ 



+ 



-A. 



2 1 — iuT]^ 1 + m?7_ 

Numerically, we make use of the cumulative function and the expected value of the variation of 
F up to some strike K. They can be computed directly by FFT techniques, without integrating the 
distribution. However one has to introduce an auxiliary distribution with known Fourier transform, 
such as a Gaussian distribution and compute corrections to this distribution. 

A cumulative function has limit in — oo and 1 in +oo: its Fourier transform is not well-defined as 
a function. Subtracting the cumulative distribution of a gaussian distribution, the function which has 
to be transformed has limit zero in both — oo and +00 and can be considered to live on a circle. FFT 
gives the Fourier transform of this difference of cumulatives which has to be added to the cumulative 
function of the gaussian distribution to recover the cumulative distribution of the jump process. In 
fact, the difference function obtained by FFT is obtained on a grid; a cubic spline interpolation is 
used to get values everywhere. We take a gaussian with same mean and variance as the process 
considered in order to reduce numerical errors. 

The second function we make use of in the computation is the expected value of the variation of 



-F up to some strike K, E 



It may be rephrased as a cumulative function on a Kou 



process with a different set of parameters and it can thus be obtained by the method just described. 
The new set of parameters is 



/ 7 \ 

a 

A± 
\r]± J 



( i \ 

a' 
A'± 



/ 7 + ^' \ 
a 

A±(lTr/±)-i 



In the case where both \± vanish, ri± are no longer relevant, this reduces to the usual Black-Scholes 
computation. 



Appendix B Closed formulas 

In the following, we consider a time-homogeneous vanilla CPPI with equally spaced rebalancing dates: 
<t\ < ■ ■ ■ <tn with li — ti-i = T for i = 1, . . . , n. From equation ([3| and from the definition of 
the standard exposure ([2]), we obtain 



Xk+i = m{Xk - ir ^ - m{Xk - 1)^ + Xk 



or, equivalently 



{Xk - 1) 



m%ti _ _|_ 1 



if ^fc > 1 



Xk-l, if^fc<l 



This formulation is useful when seeking information on the sign of Xk — 1 which enables us to tell 
wether the CPPI fell below the threshold or not. Indeed if we denote by S the (stochastic) index 
defined as 

S = min {A; G {1 . . . n} such that X^ < 1} 
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the previous relation could be rewritten as 



SAfc+l 



Xfc+i - 1 = - 1) n 



1=1 



F, 



m- 



F, 



m + 1 



(8) 



We define the local short fall probability as P^^^ = P ^Xi < 1 | Xq > l] (recall that the problem 
is invariant by translation in time), we have 



5LSF 



Fi m — 1 

— < 

Fq m 



The gap proportion (or short fall probability) is given in term of the local short fall probability as 

= P [X„ < 1 I Xo > 1] = 1 - (1 - PLS^)" 
The (undiscounted) price of a European put option with strike G is given by 
Put(to) = GE[(1-X„)lx„<i] 

= G(e [(1 - Xn)] - E [(1 - Xn) lx„<l] ) = - Xo - E [(1 - Xn) ls=oo] 

From equation ([8|), we have 

F. 



E [(1 -X„) 15=00] = (l-Xo)E 



15=00 n 



i=l 



m- 



m + 1 



where we can explicitly decompose the event {S = 00} as 



i=l 



The increments being stationary and independent, we get 

E [(1 - X„) 15=00] = (l-^o)^" 

with 

Fi 



A 
B 



E 
E 



m— - m + 1 1fi ^-i 

-TO J Fq^ m 

I<0 J Fq- ra 



1 - B 



mE 



Fx 



1fi 



"0 fQ- "1 



+ (1 - m)P 



LSF 



The final result for the put reads 



Put(to) = G(l-Xo)(l-^") 



When the CPPI is already started with to ^ ^ < ^ij separating the contribution from the first period 
yields the following formula 

Put(t) = GE [(1 - Xi) lxi<i] + GE [(1 - Xi) lxi>i] (1 - 

Using equation (|3]), one gets 

^Fx 



[Xi < 1 1 Xo > 1] 



Ft 



< K 
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and 



where 



^0 



Fi 



1fi 



<K 



K 



Ft Ft 

FqITI — \ 

Ff m 



+ (1 - Wo)Xq 



Fi 



< K 



In the case of a Black-Scholes model, one gets from appendix 



E 



that 



M{d+) and 



AA((i_). The undiscounted put price is then exphcitly given by 

Ft 



Put(t) = G(AA(d+) - AA(d_))^"-i + g{i - w^Xoy - (1 - w^)X^) (l - A^-^) 



To ihustrate the behaviour of the CPPI gap risk, we plot in figure 11 the BS price of a vanilla CPPI 



put (the value of the guarantee) depending on the risky asset spot around the current spot value. 
Note the sign change in the Delta, with high Gamma. 
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Figure 11: Price of a CPPI Put depending on the underlying spot, for the vanilla CPPI described in 
section^ The curve is asymptotically affine in both directions, with a much higher slope on the left, 
where the Put is in the money. 
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